Solitons, peakons and periodic cusp wave 
solutions for the Fornb erg- Whit ham equation 



Jiangbo Zhou*, Lixin Tian 

Nonlinear Scientific Research Center, Faculty of Science, Jiangsu University, 
Zhenjiang, Jiangsu 212013, China 



Abstract 

In this paper, we employ the bifurcation method of dynamical systems to investigate 
the exact travelling wave solutions for the Fornberg-Whitham equation ut — u xx t + 
u x +uu x = uu xxx +3u x u xx . The implicit expression for solitons is given. The explicit 
expressions for peakons and periodic cusp wave solutions are also obtained. Further, 
we show that the limits of soliton solutions and periodic cusp wave solutions are 
peakons. 
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1 Introduction 

The Fornberg-Whitham equation 

u t ~ u xxt + u x + uu x = uu xxx + 3u x ll xx , (1-1; 



has appeared in the study of qualitative behaviors of wave breaking [1,2]. It 
is a nonlinear dispersive wave equation. Since Eq.(l.l) was derived, little at- 
tention has been paid to studying it. In [3], Fornberg and Whitham obtained 



a peaked solution of the form u(x,t) = Aexp(— | 
arbitrary constant. In [4], we constructed a type of 
solutions for Eq.(l.l), which are called kink-like and antikink-like wave so- 
lutions. Unfortunately, the results in [3,4] are not complete. In the present 
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paper, we continue to derive more travelling wave solutions for Eq.(l.l), so 
that we can supplement the results of [3,4]. 



The remainder of the paper is organized as follows. In Section 2, we discuss the 
bifurcation curves and phase portraits of travelling wave system. In Section 3, 
we obtain the implicit expression for solitons and the explicit expressions for 
peakons and periodic cusp wave solutions. At the same time, we show that 
the limits of solitons and periodic cusp wave solutions are peakons. A short 
conclusion is given in Section 4. 



2 Bifurcation and phase portraits of travelling wave system 

Let u = with £ = x — ct be the solution for Eq.(l.l); then it follows that 

— ctp + cp + if + (ftp = pp + Sp p . (2.1) 

Integrating Eq. (2.1) once we have 

p"(p -c)=g-cp + p + ^p 2 - (p')\ (2.2) 

where g is the integral constant. 

Let y = p'; then we get the following planar dynamical system: 



dp 

dy g-ap + p+\p 2 - y 2 
d\ p-c 



(2-3) 



with a first integral 

H(p, y) = (p- c) 2 (y 2 - \p 2 + { l -c - 2 -)p + -l c -g) = h, (2.4) 



where h is a constant. 

Note that (2.3) has a singular line p = c. To avoid the line temporarily we 
make transformation d£ = (tp — c)d(. Under this transformation, Eq.(2.3) 
becomes 

(p - c)y 



(2.5) 

dy ,,122 
= g - C p> + p + 2<p z - y. 
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System (2.3) and system (2.5) have the same first integral as (2.4). Conse- 
quently, system (2.5) has the same topological phase portraits as system (2.3) 
except for the straight line <p = c. Obviously, ip = c is an invariant straight-line 
solution for system (2.5). 

For a fixed h, (2.4) determines a set of invariant curves of system (2.5). As h 
is varied, (2.4) determines different families of orbits of system (2.5) having 
different dynamical behaviors. Let M(ip e ,y e ) be the coefficient matrix of the 
linearized version of (2.5) at the equilibrium point (ip e ,y e ); then 



M(p e ,y e )= ye ^ (2.6) 




and at this equilibrium point, we have 

J(<p e , y e ) = det M(ip e , y e ) = -2y\ - (tp e - c)[p e - (c - 1)], (2.7) 
p(<p e , y e ) = trace(M(v? e , y e )) = -y e - (2.8) 



By the theory of planar dynamical systems (see [5]), for an equilibrium point 
of a planar dynamical system, if J < 0, then this equilibrium point is a saddle 
point; it is a center point if J > and p = 0; if J = and the Poincare index 
of the equilibrium point is 0, then it is a cusp. 

By using the first-integral value and properties of equilibrium points, we obtain 
the bifurcation curves as follows: 

gi(c) = \(c-iy, (2.9) 
92(c) = i( c - 1)*--L, (2.10) 
93(c) = i(c - l) 2 - \. (2.11) 



Obviously, the three curves have no intersection point and g^(c) < (72(c) < 
<7i(c) for arbitrary constant c. 

Using the bifurcation method for vector fields (e.g., [5]), we have the following 
result which describes the locations and properties of the singular points of 
system (2.5). 

Theorem 2.1 For given any constant wave speed c ^ 0, let 

Pi± = c-l±^(c-iy-2g for g < 9l (c), (2.12) 
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y 1± = ±Jg - -c 2 + c for g>g 3 (c). 



(2.13) 



Then we have 

(1) If g < (73(c), then system (2.5) has two equilibrium points (<^i_,0) and 
(<Pi + ,0), which are saddle points. 

(2) If g = g-s(c), then system (2.5) has two equilibrium points (c — 2,0) and 
(c, 0) . (c — 2, 0) is a saddle point and (c, 0) is a cusp. 

(3) If gz(c) < g < 92(c), then system (2.5) has four equilibrium points (<£>i_,0), 
(y?i+,0), (c, ?/!_) and(c,yi + ). (</?]_, 0) is a saddle point and (ipi + ,0) is a center 
point enclosing the orbit which connects the saddle points (c, yi~) and (c, yi+). 

(4) If 9 — 92(c), then system (2.5) has four equilibrium points (c — §,0), (c — 
|, 0), (c, — |) and (c, |), which satisfy H(c — |, 0) = H(c, — |) = H(c, |) and 
form a triangular orbit which encloses the center point (c — 3, 0). 

(5) If 92(c) < g < gi(c), then system (2.5) has four equilibrium points (y?i_,0), 
(</?!+, 0), (c, yi_) and (c, (<^ 1+ , 0) a center point enclosing the orbit 
which is homoclinic for the saddle point ((pi-, 0). 

(6) If g = Qi(c), then system (2.5) has three equilibrium points (c — 1,0), 
(c, —^y) and (c, ^). (c — 1,0) is a cusp, (c, — ^) one? (c, ^) are iwo saddle 
points. 

(7) If g > <7i(c), inen system (2.5) has two equilibrium points (c, and 
(c, . Taey are saddle points. 

The phase portraits of system (2.5) are given in Fig.l. 



3 Solitons, peakons and periodic cusp wave solutions 

Suppose that = x — ct) is a travelling wave solution for Eq.(l.l) for 

£ G (—00, +00), and lim </?(£) = A, lim </?(£) = B, where A and I? are 

two constants. If A = B, then <p(£) is called a soliton solution. If A ^ B, 
then </?(£) is called a kink (or an antikink) solution. Usually, a soliton solution 
for Eq.(l.l) corresponds to a homoclinic orbit of system (2.3) and a periodic 
travelling wave solution for Eq.(l.l) corresponds to a periodic orbit of system 
(2.3). Similarly, a kink (or an antikink) wave solution of Eq.(l.l) corresponds 
to a heteroclinic orbit (or the so-called connecting orbit) of system (2.3). The 
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(a) (b) (c) (d) 




(e) (f) (g) 

Fig. 1. The phase portraits of system (2.5). (a) g < 53(c); (b) 5 = 33(0); (c) 
53(c) < g < c/ 2 (c); (d) g = 52(c); (e) 52(c) < 5 < 51(c); (f) 5 = 51(c); (g) 5 > 51(c). 

graphs of the homoclinic orbit, periodic orbit and their limit cure are shown 
in Fig.2. 

The following lemma gives the relationship of soliton solutions of Eq.(l.l) and 
homoclinic orbits of system (2.3). 

Lemma 3.1 Assume that T is a homoclinic orbit of system (2.3) and its 
parameter expression is tp = and y = y(£)/ then u = with £ = x — ct 
is a soliton solution for Eq.(l.l). 

Proof. From Fig. 1(e), we can see that the homoclinic orbit V encloses 0) 
and connects (</?i_,0). Therefore, lim^i^oo ip(£) = ifi_. 

On the other hand, u = </?(£) is the solution for system (2.3). This implies 
that u = </?(£) is the solution for Eq.(2.1). Thus, u = (p(x — ct) is the soliton 
solution for Eq.(l.l). ■ 




(a) (b) (c) 

Fig. 2. The orbits of system (2.3). (a) The homoclinic orbit (corresponding to 
52(c) < 5 < 51(c)). (b) The limit curve of homoclinic orbit and periodic orbit (corre- 
sponding to 5 = 52(c)). (c) The periodic orbit (corresponding to 53(c) < 5 < 52(c)). 
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In Fig. 2(a), the homoclinic orbit of system (2.3) can be expressed as 
,(<p- ^i_)vV 2 + hp + k 



2(V> - c) 

where 



for ip!_ <cp< ip 2+ , (3.1) 



h = ~(1 - 3c - 3yJ(c-iy-2g), (3.2) 



l 2 = ~(1 - 4c + 3c 2 - 3#+ (3c + l)y/(c - l) 2 - 2^), (3.3) 



tp 2+ = ~{l -3c- 3y/(c-iy-2g + 2 



Substituting Eq.(3.1) into the first equation of system (2.3) and integrating 
along the homoclinic orbits, we have 

ds = --{£]. (3.5) 



■p 



[s - (p 1 -)\ / s 2 + hs + l 2 2 

It follows from (3.5) that 

/% 2+ )=/%)exp(-^|), (3.6) 



where 



W (2y/aiy/<p* + ^ + Z 2 + hp + / 3 ) ai ' 1 ' ' 



l 1 = ~(l-3c-3^(c-iy-2g), (3.8) 



Z 2 = -(1 - 4c + 3c 2 - 3g + (3c + l^c - l) 2 - 2g), (3.9) 



/ 3 = -(2 - 5c + 3c 2 -6^+ (3c + 2) v /(c - l) 2 - 2g), (3.10) 



ai = 4(1 - 2c + c 2 - 2g + J(c - l) 2 - 2g), (3.11) 



&i = -3 -4 v /(c-l) 2 -2( 7 , (3.12) 



QI = _ ; 1+v(c - 1)2 - 2g _ 



2,/(c-l) 2 -2 9+v /(c-l) 2 -2 9 



(3.6) is the implicit expression for solitons for Eq.(l.l). We show the graphs 
of the solitons in Fig. 3 under some parameter conditions. From Fig. 3, we can 
see that when 52(c) < 5 < 51(c) and g tends to 52(c), the solitons lose their 
smoothness and tend to peakons. 




Note the following facts: when 52(c) < 5 < 51(c) and 5 tends to 52(c), the limit 
curve of such homoclinic orbit of system (2.3) is a triangle with the following 
three line segments (see Fig. 2(b)): 

1 4 

2/ = ±2^~ C+ 3) f° r ^1- < ^2, (3.14) 



and 



for 



2 2 
- < y < -. 

3 ~ y ~ 3 



(3.15) 



Let us have 52(c) < 5 < 51(c) and 5 tends to 52(c); then we obtain that 

p(0 = 5«p(-i|£|) + C -£, (3.16) 
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which implies that for arbitrary constant c 7^ 0, Eq.(l.l) has peakons 



4 1 4 

u(x,t) = - exp(--\x - ct\) + (c- - 



(3.17) 



Obviously, u has peaks at x — ct = 0. We show graphs of the peakons in Fig. 4 
under some parameter conditions. 
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(a) (b) 
Fig. 4. The peakons for Eq.(l.l). (a) c = -1 ; (b) c = 2. 

Remark 3.1 (1) If we take c = | in (3.17), then we can see that (3.17) agrees 
with the result in [3]. 

(2) In the phase portaits, the triangle curve corresponds to a peakon solution. 

We have the following lemma, similar to Lemma 3.1, which indicates the rela- 
tionship of periodic wave solutions for Eq.(l.l) and periodic orbits of system 
(2.3). 

Lemma 3.2 Assume that T is a periodic orbit of system (2.3) and that its 
parameter expression is tp = y?(£) and y = y(£,); then u = <p(£) with £ = x — ct 
is a periodic wave solution for Eq.(l.l). 

In Fig.2(c), the periodic orbit can be expressed as 



1 _ .1 2. 1 _ 1 „ . 

2/ = ± V4^ ^2 C ~3^~4 C+ 3 C + 9 f° r ^-<^<c, (3.18) 



and 



(p = c for < y < y 1+ , 



(3.19) 



where 

f2- 



-(-4 + 30+^/2(9^-180 
o 



8-185-)). 



(3.20) 
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Substituting (3.18) into the first equation of system (2.3) and integrating alonj 
the periodic orbit, we have 



c 

I 



ds = for £ < 0, (3.21) 



y-(2c-§V-c 2 + fc + 4£ 



and 



c 

/ 



lip' 



(2c-|)^-c 2 + |c + 4^ 



-£ for £>0. 



(3.22) 



It follows from (3.21) and (3.22) that 

^(£) = Z + exp(-^|£|) + Le X p(^|£|) + (c-~) for tp 2 _ < y < c, (3.23) 



where 



?± = -(4±3^ + 4c-2c 2 ). (3.24) 



Let 



!T = 2|ln(y> 2 _-c+-)-]n(2Z_)|. 



(3.25) 



Then 



u(x, t) = ip(x -ct- 2nT) for (2n - 1)T < x - ct < (2n + 1)T, (3.26) 



are periodic cusp wave solutions for Eq.(l.l) with 2T period. Clearly, when 

93(c) < g < g 2 (c) and g -> # 2 (c), T -> 00, /+ -> |, L -> 0, and u(z,£) in 
(3.26) tends to 

u(z,t) = ^exp(-^|x - ct\) + (c - -). (3.27) 



(3.27) is identical with (3.17). The graphs of some periodic waves for Eq.(l.l) 
are shown in Fig. 5 under some parameter conditions. From Fig. 5 we can see 
that (73(c) < g < (72(c) and g tends to (72(c), the periodic cusp wave solutions 
tend to peakons. 



9 



-6 -4 -2 2 4 6 



-10 -5 5 10 



(a) 



(b) 



1.8- 




1.6 




1.4 




1.2 




1.0 













1.8- 




1.6- 




1.4- 




1.2 








{ i 


I i 



-80 -60 -40 -20 20 40 60 80 



(c) (d) 

Fig. 5. The periodic cusp wave solutions for Eq.(l.l). (a) c = 2, g = 0.3; (b) c = 2, 
g = 0.4; (c) c = 2, g = 0.4444444 ; (d) c = 2, g = 0.444444444444. 

Remark 3.2 (1) In the phase portaits, the semi- elliptic closed curve with one 
side on the singular line tp = c corresponds to a periodic cusp wave solution. 

(2) In [4], we dealt with the case g < gi(c), and obtained the kink-like and 
antikink-like wave solutions for Eq.(l.l). 



4 Conclusion 



In this work, by using the bifurcation method, we obtain the analytic ex- 
pressions for solitons, peakons and periodic wave solutions for the Fornberg- 
Whitham equation, given as (3.6), (3.17) and (3.26), respectively. We also 
show the relationships among the solitons, peakons and periodic cusp wave 
solutions. 
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